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Abstract 

When some granular material contained into a silo is pushed upwards with 
a piston, an irregular stick-slip motion of the system of grains is observed. We 
show how one can adapt the 'Scalar Arching Model' (sam) - proposed as a 
model for giant stress fluctuations in silos - in order to describe this stick-slip 
phenomenon. As a function of the sensitivity of the system to mechanical noise, 
the system exhibit two different phases: a 'jammed' phase, and a 'sliding' phase 
where irregular stick-slip is observed. We analyze the transition, which is found 
to be of mean-field type, and study the statistical properties of the intermittent 
stick-slip motion. 



Stick-slip motion is a very common phenomenon which occurs when two solids 
slide on each other, and has been much studied in the recent years in connection with 
solid friction [|l], Stick-slip motion also occurs in granular materials [|3], f|, ||. For 
example, the Jussieu group has performed the following experiment 0, a 2D vertical 
cell containing aluminium beads is pushed upwards by a piston through to a spring 
tightened at a constant speed. A very irregular stick-slip motion is observed, which 
can be attributed to the fact that force propagation in granular media is strongly 
inhomogeneous: most of the force is concentrated on stress paths (or arches). Some 
configurations of these paths propagate all the extenal pushing force F to the walls, 
and jamming occurs. Other configurations propagate the pushing force right up to 
the free surface, and the beads move upwards (slip). The temporal fluctuations of 
these stress paths then lead to an intermittent stick-slip motion, which we attempt to 
model here using the 'Scalar Arching Model' (sam). This model was proposed to model 
giant stress fluctuations in silos The SAM is an extension of the Chicago group's 
stochastic model |7|, which allows for the formation of arches. We have shown how 
small perturbations could lead to sudden transformations of these arches, and thus to 
large fluctuations of the weight on the bottom plate of a silo. The idea is to adapt the 
same model to the experimental situation described above, and to quantify how the 
formation of arches can indeed lead to very irregular stick-slip motions. 
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Figure 1: Beads are confined into a 2D silo. An upwards force F is applied on a piston at the 
bottom of the silo. Such a system of beads gives rise to an irregular stick-slip motion. In the 
SAM model, beads propagate forces according to rules which combines random propagation 
(encoded by the random numbers q±) and arch formation. In all the present paper, simulations 
have been performed with L = 30 (bead radii), which is close to the experimental system of 



Let us briefly recall the main features of the SAM. The granular packing is repre- 
sented in 2D by a regular lattice: each site is a 'grain' labelled by two integers (i,n) 
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giving its horizontal and vertical coordinate. We neglect the weight of the grains com- 
pared with the applied force F, and focus on the transmission of F through the grains 
from the bottom of the cell (n = 0) to the walls (i = ±L) and the free surface (n = H), 
see figure [3]. Each grain then supports the force w of its two downstairs neighbours, and 
shares its own load randomly between its two upstairs neighbours. The corresponding 
scalar equation for force propagation is thus: 

w(i, n) — q + (i + 1, n — l)w(i + 1, n — 1) + g_(i — 1, n — l)w(i — 1, n — 1) (1) 

q±(i,n) is the fraction of the force transmitted to the the grain (i =p l,n + 1), and is 
a random variable between and 1, subject to the conservation constraint q + (i,n) + 
q^(i,n) = 1. These random coefficients model the randomness of the local packing, size 
and shape of the grains, etc. At this stage, the model is the one considered in ||, 
albeit upside down. We now include a 'local slip condition': when the shear on a given 
grain is too strong, the grain can slip and lose its contact with its neighbours opposite 
to the direction of the shear. More precisely, we introduced a threshold R c such that 

g + (i,n) = l-g_(i,n) = if W ~ ~ W . + > R c (2) 

w(i, n) 

q-(i,n) = 1 -q+(i,n) = if W+ . W ~ > R c (3) 

w(i, n) 

where w± = q±(i ± 1, n — l)w(i ± 1, n — 1). These rules lead to arch formation 0. 

The walls play a crucial role in the stick-slip process, since the friction forces there 
can balance the external force F and allow the system to jam. We thus introduced a new 
parameter, the 'jamming ability' a, such that with probability a the load w(±L,n) is 
completely 'absorbed' by the wall which balances all the force carried by the grain. With 
probability 1 — a, we apply the same rule as in the bulk, i.e. for i = ±L, the fraction 
q T (±L,n) of the load w(±L,ri) hits the wall. In addition to the random q±(i, n), we 
thus define 2(H + 1) uniform random numbers at±(n) which, compared to a, decide 
whether the site (i = ±L,n) is 'absorbing' (i.e. if the grain is supported only by the 
wall) or not. In the following, the g±(z,n)'s are also chosen uniformly between and 1, 
although other choices lead to the same qualitative conclusions. 

In order to reproduce the overall stick-slip motion of the assembly of grains, we 
propose to capture the different physical phenomena which occur as follows. For a 
given applied external force F, and a given set of random numbers q±(i, n) and a±(n), 

• we calculate the total forces on the walls F w and on the free surface of the silo Ff s . 
Obviously, F = F w + F fs . 

• if Ff s = 0, the grains do not move, corresponding to a stick situation. We then increase 
the applied force F by some fixed amount AF and the time t by At. In order to mimic 
the mechanical noise which necessarily occurs when the external load is increased and 
might trigger some local rearrangements, we also change a fraction p of all the random 
numbers and recalculate Ff s . 

• if Ff s > 0, the equilibrium condition for the top grains is not satisfied, which means 
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that grains are moving. It is a slip situation. Correspondingly, the spring loosens and 
the applied force F is decreased by AF. We also change all random numbers (because 
the flow motion completely rearranges the packing). The flow stops when a randomly 
chosen configuration has enough 'anchoring' sites at the walls to yield Ff s = 0. 

The simulation starts at t — with F = and lets F increase progressively. It is 
important to note that our model is actually purely static: no dynamics is explicitly 
included. Therefore, the motion of the grains during a slipping event is assumed to be 
infinitely quick on the scale of sticking events (t is thus kept constant during slipping 
events). We thus actually describe only sticking situations, separated by slipping events 
which have two effects: untighten the spring governing the external force F, and reini- 
tialize the structure of the packing (i.e. the random numbers). As seen on figure || such 
an 'algorithm' indeed leads to an irregular stick-slip motion. Note that for a = the 
probability that Ff s vanishes is exponentially small in H. Physically, this means that, 
in order to resist to the external force F, the system of beads must generate arches 
which are strongly 'anchored' by the walls. 

The model is controlled by four parameters. The first two - namely R c , the threshold 
of the SAM, and b the aspect ratio of the cell - are of secondary importance: they do 
not affect the general features of our results. On the other hand, the jamming ability 
a of the walls and the role of the mechanical noise which modifies the local structure 
of the packing, measured by p. For R c and b fixed, depending on the values of a and p, 
two distinct phases are found. For small p or large a, the system is jammed in the sense 
that although the grains move from time to time, F goes to infinity as time increases. 
This phase can be described by the average rate of increase of F, s — F(t)/t. On the 
other hand, for large p or small a, the system slips very easily. This sliding phase is 
characterized by the fact that F(t) always goes back to 0. The relevant quantity in this 
phase is the delay r between two consecutive times where F vanishes. Figure |2] shows 
some typical plots of F(t) in these two phases. The transition between the two phases 
occurs for a critical value of a = a c (p), which allows us to obtain the phase boundaries, 
as plotted on figure |3|. Experimentally, these two regimes could be reached by preparing 
the system in different ways. A compact system would provide a rigid structure (small 
p). By contrast, a loose packing would be subject to large rearrangements (large p). 
Similarly, the state of the walls allows one to change the value of a. An experimental 
recordings of F(t) actually show parts of the two different regimes ||, |§. 

We studied how the system behaves near criticality. On figure we show the 
integrated histogram of r for R c , b and p fixed, and for different values of a. We see 
that r tends to be power-law distributed as a — > a c , with an exponent —1/2. Note 
that, as argued below, this power-law corresponds to the first return probability of 
a one dimensional random walk, (r) diverges for a = a c \ we found numerically that 
(r) oc l/(a c — a) for a < a c . In the same way, we found that the average slope of F 
versus time behaves like (s) oc (a — a c ) for a > a c (see figure ||). 

These critical laws can be understood within a simple mean field analysis. Neglecting 
the correlations, the temporal evolution of F can be approximated as a Markovian two- 
state process. Suppose the system is sliding at time t. We call p s the probability that 
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Figure 2: These plots show the temporal evolution of the applied force F for (a) a = 0.83 < a c 
(sliding phase) and (b) a = 0.85 > a c (jammed phase). Both plots have been obtained with 
R c = 0.5, p = 0.01 and 6 = 1 for wich a c ~ 0.838. 
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Figure 3: The curve a = a c (p) separates the sliding phase (below) from the jammed phase 
(above). The phase diagram (a) has been plotted for R c = 0.5 and two aspect ratios, 6 = 1 
and 6 = 3. Alternatively, we can represent the phase diagram by the curve a = a c {R c ). The 
parameters chosen for the figure (b) are p = 0.1 and again, 6=1 and 6 = 3. 
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Figure 4: These curves represent integrated histograms of the first return time r, i.e. the 
interval of time between two times where F vanishes. They have been computed with R c = 
0.5, p = 0.01, 6 = 1, and with differents values of a indicated on the plot. As a — > a c , this 
histogram gets broader and broader, and tends to the power law r -1 / 2 which is characteristic 
of the return time of simple random walks. 
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Figure 5: Below the transition, the system is characterized by the averaged first return time 
(r), and above it by the averaged slope (s) of the applied force F versus time. Near the 
transition, we find that (r) diverges like l/(a c — a) and that (s) grows like a — a c . This plot 
has been computed with R c = 0.5, p = 0.01 and 6 = 1. Linear regressions for 1/ (r) and (s) 
give respectively a c = 0.839 ± 0.004 and a c = 0.837 ± 0.004. 



it is still sliding at time t + At. Similarly, we call q s the probability of sliding at time 
t + At knowing that the system is in a jamming configuration at time t. Obviously, p s 
depends on a and q s onp and a. For example, one has q s (p = 0) = and q s {p = 1) = p s . 
Figure^ shows q s {p), as determined numerically. This simple two-state model can be 
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Figure 6: The probability q s that a system in a jamming configuration at time t becomes 
sliding at time t + At, depends on p. For fixed R c , b and a, we call p c the critical value 
of p such that a = a c . Three points of the curve q s {p) are analytically known: q s (0) = 0, 
q s (l) = p s and q s (p c ) = 1 — Ps- The whole curve has been obtained numerically for R c = 0.5, 
a = a c (p = 0.1) = 0.914 and 6=1. 



explicitly solved. The critical point is found to be when the probability of sliding after 
jamming is equal to the probability of jamming after sliding, i.e. q s (p) = l—p s . At this 
point, the probability that F increases is equal to the probability that F decreases, 
which implies indeed that F behaves as a random walk. Provided that the functions 
p s and q s are regular near the critical line, one also finds the observed linear behaviour 
of s and 1/ (r). The fact that the temporal correlations are found to be small however 
means that our model fails to capture 'precursor' effects before the slip, which have 
been observed experimentally in || . This is related to our simple rule where each grain 
can 'move' under the influence of the external noise with equal probability. Finally, 
we also looked at the two following quantities: the distribution of the heights of the 
slips, and the distribution of the intervals of time between two slips. The tails of these 
distributions are found to be exponentially decaying. 

In conclusion, we discussed how the Scalar Arching Model can be modified to de- 
scribe granular systems undergoing intermittent stick-slip motion. We showed that 
such a system can present two different phases: a 'slipping' stick-slip (or sliding) phase, 
where the external pushing force remains finite, and a 'sticking' stick-slip (or jammed) 
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phase. The transition between these two regimes is found to be of mean-field type. 
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